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We calculate the deflection of a light ray caused by the gravitational fleld of a cosmic string loop 
in the weak field limit and reduce the problem to a single quadrature over a time slice of the loop's 
world sheet. We then apply this formalism to the problem of gravitational lensing by cosmic string 
loops. In particular, we find an analytic solution for the special case of a circular loop perpendicular 
to the optical axis. As examples of more complicated loops, we consider two loops with higher 
' frequency Fourier modes. The numerical analysis illustrates the general features of loop lenses. Our 

I estimates, using typical parameters for GUT scale loops, show that the stringy nature of loop lenses 

^\ . can be observed for lensing systems involving high redshift galaxies {z ~ 2), and we suggest that 

gravitational lensing can confirm the existence of GUT scale strings if they are the seeds for large 
scale structure formation. 
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I. INTRODUCTION 



Current research focuses on two scenarios for the formation of structure in the universe: the first where structure 
formation was seeded by adiabatic perturbations produced during an inflationary epoch and the second where structure 
accretes around isocurvature perturbations produced by topological defects such as cosmic strings, global monopoles 
or textures. In the latter scenario it should be possible to directly detect the presence of topological defects in the 
■ present universe leading to immediate confirmation of the scenario. On the other hand, the lack of direct evidence for 
^3 [ topological defects in the present universe can lead to constraints on the defect scenario for structure formation and 
0^ perhaps be considered as circumstantial evidence in favor of the inflationary alternative. Thus it is quite important 
[ to consider specific distinctive signatures of the various topological defects that can be used to directly observe them. 
Let us specifically consider cosmic strings, the model which will be relevant to the work in this paper (for a 
review of cosmic strings, see Ref. [Q). A number of observable features produced by cosmic strings of mass density 
suitable for structure formation have been discussed in the literature. These include discontinuous patterns in the 
microwave background radiation 0], generation of a gravitational wave background that could be detected by noise 
in the millisecond pulsar timing ^ and gravitational lensing |0,^,^;0J|]- The ongoing observations of anisotropics 
in the microwave background radiation are expected to yield stronger constraints or positive results over the next 
decade or so. The millisecond pulsar observations can only impose tighter constraints on the string scenario since a 
positive detection of gravitational waves does not specifically imply the existence of cosmic strings. There has also 
] been sporadic effort over the last decade to work out the gravitational lensing signature of cosmic strings but, perhaps 
due to the difficulties encountered in understanding the evolution of the string network, no distinctive result emerged 
from these analyses. However, an analytical framework for describing the string network has been constructed over 
the last few years and the time seems ripe to reconsider gravitational lensing as a tool for searching for strings. The 
timing is also right from the observational viewpoint since several new initiatives are underway that promise to survey 
much wider and deeper regions of the sky. p. 

In this paper we investigate gravitational lensing by cosmic string loopstl. We begin in Sec. II by estimating the 
probability of string lensing and in so doing we review some of the relevant properties of cosmic strings. Our estimates 
are based on recent results for the string network evolution summarized in Ref. 0]. Next, in Sec. Ill, we consider 
photon propagation in the metric of a loop. The problem appears to be quite difficult at first because the oscillating 
loop is a complicated time dependent gravitational source. Yet we are able to show that the problem reduces to one 
that is static where a specific time slice of the loop's world sheet is suflicient to determine the gravitational lensing 



'^Gravitational lensing by global monopoles and textures is likely to be less interesting since only a few of these are expected 
to occur within our horizon. Also, their spherical symmetry will lead to lensing that is harder to differentiate from that due to 
conventional sources. 



1 



effects. In other words, the bending of hght by a string loop is equivalent to the bending of light by a static curved 
rod with non-uniform energy density, a result similar to that for the energy shift of a photon propagating in a string 
loop background first derived by Stebbins We also rederive the energy shift of the photon in the Appendix and 
recover a logarithmic term that appears to have been eliminated by the regularization procedure used in Rcf. 

Once we h ave s et up the formalism for an arbitrary loop and described some rudiments of gravitational Icnsing 
theory (Sec. IV A ) , we apply it to treat the Icnsing due to a circular loop that is oriented in a plane normal to the 
optical axis (Sec. IV B). The results for the circular loop are in agreement with the assumption in Rcf. [|| that 



the photons passing through the loop remain undeflected. The deflection of a photon trajectory not threading the 
loop can also be described quite simply and the whole problem can happily be solved by hand without resorting to 
numerical evaluation. 

The perpendicular circular loop, however, is a very special case as even a change in the orientation of the loop 
yields qualitatively different results, and loops with less symmetry have completely different Icnsing behavior. The 
assumption that photons passing through the loop remain undeflected fails for general loops. We study the lensing due 
to several generic loops numerically and provide image ma ps tha t promise to distinguish cosmic string lensed images 



from more conventional gravitational lensing events (Sec. IV C ). Here we also show that the Einstein radius of the 



string loop is comparable to the typical loop size for any value of the string tension and so the stringy nature of the loop 
plays a crucial role in determining the structure of the lensed images. Effective techniques — for example, techniques 
that replace the string loop by a point mass plus perturbations — are unlikely to yield successful approximations 
leading us to conclude that string loop lenses ought to be observationally distinct from garden variety astrophysical 
lenses. 

In Sec. we summarize and discuss our main results. We also qualitatively discuss the effects of long strings and 
describe further work to come. Finally, Sec. ^ contains some concluding remarks. 



II. LENSING PROBABILITIES WITH STRING LOOPS 



Gravitational lensing by cosmic string loops is an interesting problem only if there is a realistic chance of observing 
a loop lens. In this section, we will estimate the typical size and number density of loops and use these values to 
determine the likelihood of observing a string loop. Our arguments will be based on scaling solutions to the string 
network evolution which are indicated by numerical simulations and semi-analytic treatments (for a review of scaling 
solutions see Ref. 0]). 

The string network consists of two components — the long (or infinite) strings and the closed string loops — 
assumed to have formed during a phase transition in the very early universe. If the strings are to seed large-scale 
structure formation, they should have a linear mass density /i ^ 10^^ gms/cm ( conveniently expressed in Planck 
units as G/i ^ 2 x 10~^ where G is Newton's gravitational constant). At any epoch in the history of the universe, 
curved sections of strings will oscillate under their own force of tension; colliding and intersecting strings will undergo 
reconnections; the strings will stretch under the influence of the Hubble expansion, and the oscillating strings will 
lose energy primarily to gravitational radiation. The complicated evolution of the string network has been studied in 
a number of works and a consensus is emerging that the long strings obey a scaling solution [|io|,0. In other words, 
the energy density in the long strings scales with time as ~ and the typical distance between strings also 
scales with time, L on t, where there are on the order of 10 long strings per horizon at any epoch. The reconnections 
of long strings and large loops are found to copiously produce small loops. But there is generally less agreement on 
the precise number density of loops present at any epoch. If the typical loop produced is assumed to have a length 
£ at, then numerical results (that do not take the gravitational back reaction into account) place an upper limit of 
a < 10~^ [ [lO| , [ri| . To get a lower limit on a, one supposes that highly curved sections of string straighten out on very 
short time scales due to gravitational radiation which occurs at a rate 

where, F is a numerical factor depending on the shape of the loop and has been evaluated to be ~ 60 by considering 
several classes of loops. Since the size of the loops is related to the curvature of strings, the smallest loops that can be 
produced are those that have lifetimes longer than the Hubble time scale. This leads to a « FG/i ~ 10~^ which we 
will use for the purpose of numerical estimates. Then, for a scaling solution, the number density of loops with length 
between i and £ + d£ at a given time t is 
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where v is about 0.5 ||. Now we can estimate the solid angle of sky coverage necessary to observe a string loop lens. 
Using eq. (H), the number of loops of length between at and fiott — supposing that loops larger that jiat will not 
make interesting lenses — per solid angle at redshifts less than zq is 

m_ n , f 4 dnii,t) 1 ji + z-ymy 

where t = (2/3)iJo"^(l + z)-^/'^. The various explicit factors of the redshift arise from the changing volume element 
in an expanding universe. One can evaluate this integral analytically and find the result 

dNi 27 ua([i-l) r, /, x . , .i 

-dn-Y iPa + rG,)ia ; TG,) + + " + " 

For zq = 1, f3 ^ 2 and a = TG^, there are about 400 {TG^i)"^ loops per steradian which corresponds to about 1 
loop for every (3°)^rG/X4 of sky, where we have defined rGfj.4 — TG^jl x 10''. But, this does not necessarily tell us 
how many loop lenses we can expect to observe because we have not yet factored in the odds that a galaxy will be 
near enough to a given loop to be significantly lensed. To resolve this problem we require the luminosity distribution 
of high redshift galaxies. 

For an order of magnitude estimate of the luminosity distribution of high redshift galaxies, we can use information 
about nearby galaxies and extrapolate to higher redshifts. Such an extrapolation is justified only if the effects of 
evolution can be ignored and we shall assume that this is the case. The local luminosity function, (/>, is reasonably 



approximated by the Schechter function |12|, 



= 0,w"e-", (5) 

where u is proportional to the luminosity, a = 1.07 ± .05, — Q.QlQe^-^h^ Mpc^'^ and h is the Hubble parameter in 
units of 100km s^^ Mpc^"'^. The number density of galaxies, Ng, per unit solid angle per unit redshift per unit flux is 
given by 

d^ N 

" = inH^'^zUi'iz^), (6) 



dildzdh 



where k is proportional to the flux. For a flux limited or equivalently a magnitude limited survey, we integrate over 
K from a lower bound to infinity. The lower bound can be expressed in terms of the apparent magnitude m using the 
relation 

K = 9.0 X 10^^/1^210*=^*-™, (7) 

where — —19.53 ± .25 + 5 log h is the absolute magnitude of a characteristic galaxy and a reasonable value for the 
magnitude limit is m = 24. For a flat universe with no cosmological constant, age constraints suggest that we should 
use a value oi h = 0.5 which gives a lower bound of k = 0.057. We can express the integral over the flux in terms of 
an incomplete gamma function which gives us 

^ =2.7 X 10Vr(a + 1,0.057^2). (8) 
dildz 

Let us suppose that we have a loop located at about a redshift of one away from us and we are observing galaxies 
that are about a redshift of two to three. Then integrating eq. (||) tells us that there are about 1.4 x 10^ galaxies per 
steradian to work with. The typical angular separation between galaxies is about 5.6 arc seconds, which we shall see 



in Sec. [V C is small enough that a foreground string loop can always be expected to have a galaxy in its background. 



Therefore every loop has a good chance of being seen as a gravitational lens. 



III. GEODESIC DEFLECTION BY A STRING LOOP 



In this section, we consider the deflection of null geodesies in the presence of a cosmic string loop, which in the 
limit of geometric optics, will correspond to the photon paths. (Readers only interested in lensing applications might 
wish to proceed directly to the final result given at the end of the section.) We shall assume that the background 
space-time is flat Minkowski space, and ignore the effects of curvature or universal expansion. Since we are interested 
in loops with radii much smaller that the horizon or curvature scale, this approximation will be valid for determining 
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the photon deflection. We shall further assume that the string contribution to the metric is weak and may be treated 
as a perturbation on the flat space. Thus we may write the complete metric as 



where we use the convention rj = diag(— 1, 1, 1, 1), and h^^, is a small perturbation, so all terms of O^h?) will be 
ignored. To fix the gauge, we apply the harmonic condition, that is <?^T^j^ — 0, which leaves us with a simple wave 
equation for the metric 

□ ^ V = -167rG5^„ (10) 

where 5*^1, is related to the stress-energy tensor by S^y = T^i, — l/2ri^yT^. It should be noted that for a loop of finite 
size, the metric far from the loop is asymptotically flat. 

In the geometric limit, we can treat the trajectories of photons as the null geodesies of the metric. If A is an affine 
parameter, then the momenta are given by 

The evolution of the P" are determined by the geodesic equation 

— + r;:,p^p^ = o, (12) 

with the constraint P^P^ = 0. In the weak field limit, the Christoffel symbol, T, is 

^"^u - lv"^iKs,u + Ks.^ - h^,^s)- (13) 

To zeroth order, i.e. in the absence of a loop, the Christoffel symbol vanishes, implying the momenta are constants 
and the coordinates grow linearly with A. In particular, one can choose the afhne parameter A to be t/P^ . Also, 
since the non-zero Christoffel components are all first order, we can contract T with the zeroth order momenta, and 
further consider the Christoffel components to be functions of only the zeroth order coordinates. Let us be explicit by 
considering the dimensionless four velocity which is defined by 7^ = P^ / P^^^ where P^fJ,^ is the zeroth order energy. 



It has a zeroth and first order part which we may write as 7^ — j^^^s^ + j'^^y Rewriting eq. (12) using the appropriate 
substitutions, to first order we get 

= -{h^^a^u - ^ V,")7('o)7('o)- (14) 

Let us now consider the light ray emitted from a static source residing at a distance much greater than the size 
of the loop and traveling to a distant observer, that is, we will consider our source and observer to lie effectively at 
infinity where spacetime is flat. The components of 7^ which appear explicitly on the right hand side of eq. (^^ are 
zeroth order, and the perturbation in the metric h^i, is a function only of the zeroth order coordinates. Therefore, we 
can integrate eq. ( p^ explicitly giving 

f°° 1 

However, hfia,iyP'^ = dh^^a/dt so that the first term in the integral is just a vanishing surface term. Thus only the 
second term survives and we are left with 

1 f°° 

7q(i) = ^ / dt /ip^,a7(o)7('o)- (16) 

We now require the metric produced by a string loop if we are to calculate the geodesic defiection; in the weak 
field limit, this problem has been solved. The configuration of a string loop is described by the position of the string 
/^((T, t), where t is a time variable and ct is a parameter along the loop. The equations of motion for the string in 
Minkowski space are given by 

r - n = 0, (17) 
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with the constraints 



(18) 



Dots here refer to derivatives with respect to t while primes refer to derivatives with respect to a. It is convenient to 
choose, 

f = t (19) 
and to write the string solutions as a superposition of traveling waves with non-linear constraints |[l3[| , namely, 

a((T - t) + b((T + t) 



{{a,t) = 



(20) 



a'^ = b'^ = 1 . 



(21) 



For any closed loop of length L, solutions for f must satisfy the periodic condition i{a + L,t) = f (cr, t). In the center 
of mass frame of the loop, the functions a and b are periodic as well, but this is not true in general. The description 
of the loop as given in eq. ( pO| ) also applies to loops that have a net momentum if we use the boundary condition 
a(a — t) — a((T — t + L) = b(cr + t + L) ~ b(cr + t) = A., where A is the loop's center of mass velocity. 
The energy-momentum tensor of the string is given in terms of / by 

T'"- = ^1 [ dairr - f^DS^'^ (x - f (a, t)). 



(22) 

With T^'^one can write the solution to eq. ( |l0|) as an integral over a Greens function in retarded time, explicitly giving 



us 



J |x-x' 



(23) 



where t = i — |x — x'| is the retarded time. Using the stress-energy defined in cq. (p^), we can evaluate the spatial 
integral in eq. (E3h, 



/ip^.(x,t) = iG^i da- — 
J x 



fl-(x-f)./ 



where 



(24) 



(25) 



Let us use eq. ( |24| ) to test the validity the weak field approximation. Derivatives of are of order unity, so it 
follows we may treat F^^ as order one. The integral over a gives a term of order L, the loop length. Roughly then, 
the perturbation is small, i.e. hf_,^ <^ 1, when |x — f | ^ GfiL. The string energy density is given by < 10~^, 
so that the weak field approximation will break down only for photons which come within a small fraction of the 
loop length of striking the string itself, implying that, for most lensing scenarios, the approximation is valid. There 
is, however, an exception to this estimate in the case of cusps. Cusps are points on the loop which momentarily 
achieve the speed of light as the string oscillates, and, if the velocity of the cusp points in the direction of the photon 
trajectory, the deflection can be sing ular |l|. However, cusps only occur instantaneously, and the probability of a 
photon encountering one may be neglected. 

Eq. (^4|) might seem to be the appropriate starting point for calculating ray deflections, but we shall see that it 
is more convenient to solve this problem in Fourier rather than coordinate space. To avoid confusion, let us use the 
following conventions. 



F(x'^) 



1 



27r)4 

4 ik 



d'^ke-'^^-'' 



and consider eq. (|lO|) again. Transformed, it becomes 
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Now we choose the unperturbed photon trajectory to be 

:k — jt 

and use ft.^i/,Q = —ikah^^, to write 



(26) 



I. 



dt hf_,^,a{t,jt) 



G 



dt 



S,. 



k^kx ' 



(27) 



where all integrals are implicitly evaluated over an infinite range. The integral over time can be evaluated, and 
just transforms the exponential into a delta function 2tt5(^ ■ k — /cq)- If we now decompose the wave vector k into 
components parallel and perpendicular to 7 so that k = (A:|| , k^), where k • 7 = fc|| and k_|_ is a two dimensional vector 
perpendicular to 7 then the delta function becomes 5{k\\ — ko), which allows us to evaluate the integral over ko. What 
remains is an integral over the three dimensional vector k 



2G 
2" 



TT 



la — — I dk\\d^k± 



kl 



(28) 



where k± refers to the magnitude of k^. To make any more progress, we need to evaluate S in terms of f . From eq. 
(E3) one may infer that 



(29) 



One should note that f is measured from an origin chosen somewhere along the zeroth order photon trajectory and 
not from the loop center of mass. Later we shall find it more convenient to decompose f into a component r measured 
form the center of mass and Xq measured from the origin on the photon path, but for compactness we stay with f for 
now. Replacing fco with fc|| , we can substitute eq. ( [29| ) into eq. (|2|) to get 



2Gn 



da 



Jk^ (t-/||) -iki-fi 



dt I dk\\d k±—^Ffj_^e'"»^~ ■''"e 



(30) 



where we have decomposed f into its parallel and perpendicular components. 

At this point it is necessary to choose a. If we want to know the redshift induced by the presence of a string, we 
would solve eq. ( pO| ) for ka — > fc|| ; we will do this in the appendix to compare our results with previous work, but 
for now, we are interested in the deviation of the photon path from its zeroth order direction. Thus we consider the 
components ka k±i. In this case, one may perform the integral over fcy to produce a delta function 5{t — f\i(a,t)), 
and so integrate out the time leaving 



AGfi 



da 



F„ 



d k±e' 



ik±i 



kl 



where tQ is the solution to 



/||(io,cr) = ^0- 



(31) 



(32) 



The integration over kj^ can be evaluated by recognizing that e ^^^'^^ /k^ is proportional to the Fourier transform 
of the Greens function for the two dimensional Laplacian operator, — 7rlog(/^). The final result is suprisingly simple. 



I±^ = 8Gfi / da 



Finally, if we define the two dimensional vector a. = 7j^(t — > —00) — 7j^(< 
from its zeroth order trajectory, then from eq. ([T^)and eq. (^) we will get 



a = -AGn 
where t^ is determined by solving eq. (p^ 



da 



fl 



(33) 

00) to be the deviation of the photon 

(34) 
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IV. GRAVITATIONAL LENSING 



A. Basic Lens Theory 



In the previous section we derived an expression for the deflection of the photon momentum that passes near a 
cosmic string loop. Now we would like to use this expression to determine how an intervening string loop affects the 
images of various sources. Let us define an origin which is located at the center of mass of the loop, so the vector 
r(i7, t) will trace the loop measured from this origin. Further, we shall define a second vector Xq which points from 
the loop center of mass to the location of the photon at t = 0, so we may rewrite f as the difference r — Xq. Now let 
us define the optical axis as the line connecting the loop center of mass to the observation point, that is r|| and X|| 
both point in the direction defined by the ray connecting the loop center of mass with the observer. Since the source 
rays which are likely to be affected by the loop come in at very shallow angles, we can treat them as rays parallel to 
the optical axis for the purposes of calculating a. We shall also assume that the deflection occurs instantaneously in 
the lens plane defined as the plane normal to the optical axis located at the loop center of mass. In each case, these 
approximations are accurate to terms of order the loop size over the distance of the source or observer to the lens 
(whichever is smaller). For the loops we shall be considering, these approximations will be more that adequate, since 
typical distances are on the order of the horizon while the loop size is about a factor TGfj, smaller than that. 

In Fig. we show a schematic representation of the lensing system. A source lies in a plane perpendicular to the 
optical axis a distance Di^ from the lensing plane at the point 5'. The two dimensional vector in the source plane t] 
points from the optical axis to the source. A ray emitted from the source intersects the lensing plane at /, and the 
two dimensional vector in the lensing plane ^ points from the optical axis to /. At /, the ray is deflected by the vector 
a. To be observed, the ray must intersect the observer O at a distance Di from the lens plane on the optical axis. 
For small deflection angles, this condition is met if 

r? = ^|-Asa(0, (35) 

where Ds is the distance from the observer to the source plane. This formula, referred to as the lens equation, gives 
us the location of a source, rj, given the location of its image, ^. Usually we are interested in the inverse problem, 
namely given a source location, where are its images located. Note that the inverse of eq. ( p5| ) is not necessarily single 
valued, so it is possible that a particular source has multiple images. 

We must digress for a moment to discuss exactly what is meant by distance in the lens equation, for in an expanding 
universe, several definitions of distance are possible. For example, if the universe were seeded with sources of a known 
absolute luminosity, one could define the distance to any of these sources to be proportional to the square root of 
the absolute luminosity divided by the observed flux. However, for the case of gravitational lenses (Fig. 1), we are 
interested in the angular size of an image as it appears to a terrestrial observer, and we find that the luminosity 
distance is not the best choice. We should then like to define the angular diameter distance, (Ia, as the ratio of 
the object size and the angle that the object subtends at the location of the observer. Then, using an FRW metric 
and assuming a flat, matter dominated universe with no cosmological constant, it is easy to show that the angular 
diameter distance is fl^ 



dAizi,Z2) 



Ho{l + Z2) 



1 



(36) 



The lengths Dg, Di and Dig are all angular diameter distances given by this equationH. 

Given an extended source, not only can a lens change the location and number of images observed, it can also affect 
the magnification and shape as well. To see this, let us, for the sake of convenience, first redefine the lens equation in 
terms of dimensionless variables. The loop length L is defined in terms of the energy of the loop as E/ fi. A convenient 
length scale for us is the loop radius which we define as i? = L/2it. If we define 

^4 (37) 



^A word of caution: the real universe, especially at late times, is clumped and only homogeneous on average. A better distance 
choice may come from models like the Dyer-Roeder equation [Q, but for the present we prefer not to consider these technical 
questions. 
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then the lens equation is reduced to 



a(x)=a(0^, (39) 



y = X - a(x). (40) 



One interesting question to ask, then, is given a narrow pencil beam emitted by the source which subtends a solid angle 
duj* , what is the solid angle dco subtended by its image? This question is directly related to the issue of magnification 
as the observed flux relative to the emitted flux is just the ratio of the solid angles dio/duj* . The answer may be found 
by looking at the Jacobian matrix 



which gives us the magnification factor 



A, = (41) 



and the magnification is defined as the absolute value, \fi\. For many lensing systems, there will exist curves in x 
for which /i(x) is infinite or equivalently the determinant of the Jacobian vanishes. These curves are referred to as 
critical curves, and the corresponding map of the critical curves into the source plane are called caustics. Caustics 
play an important role in determining the number of possible images. Consider a system for which the Jacobian never 
vanishes. The lens equation is then globally invertible and therefore single valued. Thus when there are no caustics, 
there can be only one image of the source. If there exists a caustic, then the lens equation is only locally invertible 
and there may be multiple images. In fact there are many general theorems which can be proved about caustics and 
images; we shall not discuss them here, but the interested reader is directed to a detailed review of the subject of 
gravitational lensing by Schneider, Ehlers and Falco [ p^ . 

We have seen how a lens can affect the magnification of an image, and we would now like to show how a lens can 
change the shape of the observed image. We shall restrict ourselves here to a discussion of small sources so that 
we may treat the problem differentially. Consider two points on a source separated by a distance Y. From the lens 
equation we can see, to first order in the Taylor expansion, that the image displacement X will be 

X = A-^Y. (43) 

Now consider a source that is a small circle. From eq. (03) we can surmise that the image will be an ellipse with major 
and minor axis pointing along the eigenvectors of ^~^vith lengths equal to their eigenvalues. The most dramatic 
results will appear near a critical curve. Here, usually one of the eigenvalues of blows up, so that images are not 
only magnified, but stretched as well. What is observed is a stretched image, but since the flux grows in proportion to 
the area, one observes an image which appears to be a stretched version of the source but with the original brightness. 



B. An Analytic Example: The Perpendicular Circular Loop 

In this subsection, we will consider a simple example of a circular loop which lies in a plane perpendicular to the 
optical axis. This problem has the virtue of being analytically soluble, and it will also illustrate some of the features 
that migh be expected for a more general loop. Most notably, there is a discontinuity in the deflection a as light rays 
go from passing through the loop to passing outside it which noticeably influences the resulting images. 

Let us begin by choosing the z direction to point along the optical axis. For the perpendicular circular loop, the 
configuration of the string is given by 

r = i?cos(-^)[cos(^),sin(^),0], (44) 

where R is the maximal radius of the loop. From eq. (^5|), one can see that for any planar loop perpendicular to the 
optical axis, ^^^'^F^y = I. Also, there is no parallel component in r, so f\\ = 0. Thus, using eq. ( |3^ ) and eq. (p7|), we 
can write the deflection as 



8 



DisDi 
RDs 



27T 



de 



^2 _|_ J.I2 ^^^^ _ 2xr' cos 6 ' 



(45) 



where r' — r(to)/i?, x = |x| and ioi defined by eq. (|3^), is a constant. The behavior of this integral is most easily seen 
by analytical continuation on to the complex plane. If we change variables to z = pe*^, then the integral in eq. (45) is 
a contour integral along the p = 1 circle in the complex plane and hence reduces to a sum of complex residues. One 
finds that for image points such that x < r' , the residues cancel and there is no deflection, but for image points such 
that X > r', the result is non-zeroEl. We point this out because in principle the deflection for any loop can be reduced 
to finding residues of a complex contour integral. In practice, however, the problem is complicated by the need to 
solve eq. ( ^2[ ) to get to as a function of a, and there appears to be no advantage in using the calculus of residues over 
numerical integration to find a, although there may be other cases aside from the circular loop where this method 
could be usefully applied. Returning to the circular loop, one finds that q:(x) = if |x| < and 



a{x.) = —SnGfi- 



DisDi x 



RD,, 



if |x| > r' 



For rays passing outside the loop, the deflection is exactly the same as if the loop were replaced by a point mass. 
Inside the loop, the ray is undeflected, a feature which distinguishes a normal circular loop from a point source. 
Given the deflection, we will now consider the magnification, specifically looking for the caustics. The Jacobian 



matrix can be written as Ay — 6ij 



a.i j , so for the circular loop we find that when x > 



A = 



1 



C{xl - x\)lx^ 
-2Cx\Xil x"^ 



1 



~2Cx\X2l x"^ 
C{x\ - xl)/x^ 



where C = SnGfiDisDi/ RDs and x is given by its components {xi,X2)- For x < r' 
det{A) = 1. For x > r' , the determinant of A is given by 

det{A) = 1 - — . 



A is just the identity matrix and 



(46) 



One real root, namely x — VC, is proportional to the Einstein radius Re = RVC. Note that this definition is 
consistent with the familiar definition of the Einstein radius of a point mass, namely R^ — \/ 'iGMDisDi/ Dg, where 
for a string the mass is just M — 27r _R/i. So, if \/C > r' , the loop lens will have a circular critical curve with radius 
-s/C, but if y/C < r' there will be no critical curve. In the latter case, the loop would be difficult to detect because 
there will be only one image of the source which will not be significantly distorted or magnified. In the former case 
{VC > r') there are, strictly speaking, two critical curves because the discontinuity at the location of the string should 
be smoothed out on a scale given by the thickness of the string. Then the first critical curve is the usual Einstein ring 
with radius \/C and the second occurs because the determinant of A must change from a positive value (namely 1) 
inside the loop to a negative value just outside the loop (namely 1 — C^/r'"') and hence deiA = when the light ray 
passes through the stringJ. For the case of the circular loop, we can analytically invert the lens equation as well, and 
after a little algebra we find 



(47) 




(48) 



when X > r' , but x = y when x < r' . To help visualize this rather complicated picture, we show in Fig. ^ a series of 
images for a circular source as it is displaced further away from the optical axis. We have chosen a loop radius r' = 1 
with a value of C = 2.25 corresponding to an Einstein radius of 1.5. It follows that there are two critical curves with 
radii = 1 and 1.5 



^This discontinuous behavior will hold for all loops and is easily understood in terms of the motion of poles in the complex 
plane. 

*This second critical curve is of no observable interest since the string thickness is a mere 10"^" cms or so. It does have formal 
significance when considering lensing theorems. In the following section, where we give image maps, we shall ignore the effects 
of light rays that actually pass through the string. 
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C. Numerical Examples 



The perpendicular circular loop is one of the few cases which may be treated analytically, but it is not a very generic 
loop. In this subsection, we would like to consider a few examples which represent more general cases of cosmic string 
lenses. We will consider loops which are neither planar nor oriented in any special way with respect to the optical 
axis and are characterized by reasonable astrophysical parameters. 

Suppose that we have a bright galaxy or cluster of galaxies at a redshift of z ^ 2 which is lensed by a cosmic string 
loop. Will a "typical" loop produce any observationally distinctive images? The answer to this question depends 
roughly on the value of the Einstein radius of the system. For a general loop, we define the Einstein radius as 



Re 



DisDi 



(49) 



so if we were to replace the loop with a point source with mass equal to the mass of the loop, the critical curve 
would be a circle with the Einstein radius. If i?e 3> i? then the images which would be significantly distorted by the 
string loop will appear to be similar to those of an equivalent point source. If i?e ^ R, then the loop will not have a 
significant affect on any images of the source as most of the mass of the loop will lie outside the Einstein radius. If 
string loops are to produce distinctive images, we should have Re ~ R. So let us consider what may be reasonable 
parameters for a loop lens system. 

We can estimate the size of a "typical" string loop on the basis of the string network evolution described in Sec. ||. 
Loops are expected to be produced from the network of long strings with typical size R ^ TG^t{z = \)/2-k. These 
loops gradually lose their energy to gravitational radiation (eq. (|l|)) but survive for about one Hubble period. (We 
are assuming that the loops are not further fragmented significantly by self-intersections. If that is the case, we would 
need to divide the size estimate by the expected number of fragmentations and factor in the survival period in the 
estimate of the lensing probability in Sec. ||.) Then, using this value for the loop radius, taking F = 60 and locating 
a source at z = 2 and the loop at 2 = 1, the Einstein radius is given by 



Rp — R\ 



F tD. 



R 



(50) 



— which is in just the right range for distinctive string lenses. Further note that the result R,, ^ R is independent of 
the string tension fi and so the stringy nature of loop lensing is important for strings of any mass density. We therefore 
expect that the images produced by string loops will have characteristic features of stringy lenses, thus offering the 
realistic hope that strings may be observed definitively through gravitational lensing U. 

We shall consider two classes of loops to illustrate the kinds of images that one may observe with cosmic string 
loops. The first is a mixture of the fundamental and the first excited mode 



R 



sin(T_i + - sin 2cr+j + [cos cr_ + - cos 2(7+]k 



(51) 



where cr_ = {a — t)/R and <j+ = {a + t)/R. We will call this the "two-loop" for short since it is a superposition of a 
fundamental mode with a frequency two mode. The second loop configuration is a class of loops found by Turok |l^] 
which is the general solution for a loop with frequency one and three Fourier modes. This solution has the form 

R ( I 1 

r= — < (1 — a) sin (T_H — a sin 3cr_ + sin fT+ i (52) 
2 1 3 

(1 — a) cos (J- + -a cos 3ct_ + cos (/)COS (t+ 
o 

2y^a{l — a) cos (t_ + sin (/)cos (t+ k| . 

We shall call this the "three-loop" for short. In general, the coordinate system which defines {i, j, k}, i.e. the loop 
coordinates, can be rotated with respect to the optical axis, which, for convenience, we shall define to be the z axis. 



^The wiggles on long strings can probably also be regarded as loops of size R and hence the lensing by long strings would 
mimick that of a linear array of loops. This would be another characteristic feature of lensing by cosmic strings. 
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The relative orientation of the systems can be described by the three Euler angles. However, we are not particularly 
interested in the orientation in the sky, so it is sufficient to consider only two of the three. Thus, if we start with 
the systems {x, y, z} and {i, j, k} aligned, we shall first rotate the loop in the x — y plane from x to y through an 
angle 9i . Then we shall tilt the loop with respect to the observer by rotating in the y — z plane from y to z through 
an angle O2, giving us the final orientation of the loop coordinates with respect to the optical coordinates. 

To find the images for a given source is a numerically non-trivial problem because one has to invert eq. (35) to 
get ^ as a function of rj. This involves finding simultaneous roots of two equations in two dimensions. But, as there 
are no general methods to solve such a problem, we are left with no choice but to map the entire image plane (^) 
into the source plane [rf) and then consider the inverse mapping. (In practice, a lattice of points in the image plane 
is mapped on to the source plane.) To locate the images for a particular source point, we test every triangle in the 
image plane that can be formed by three neighboring points on the lattice. If, when mapped to the source plane, 
the image triangle encompasses a source point, then we know that the source has an image located somewhere inside 
the image triangle. Calculating a. entails performing one numerical integration, so the inversion process is reasonably 
cheap with regards to cpu time and good resolution is possible. The triangle search method, however, has one failing 
when finding images with strings. Since a will be discontinuous for light rays which pass just to either side of the 
string, image triangles which are formed by points on either side of the loop may encompass sources which they would 
not if all the image points were on one side or the other of the loop. These are spurious images as they require the 
light rays to pass through the stri ng its elf, so we have rejected these points when constructing the images of specific 
sources (see the discussion in Sec. IV B| ). |-| 

Let us now use this technology to generate some images for the loops which we have discussed in this sectionQ. In 
Fig. H we consider the two-loop aligned with the optical coordinates so = and 02 = 90. The phase of the loop is 
given by "0 which is defined as the time at which the light ray intersects the lens plane — the plane normal to 7 and 
containing the loop's center of mass. For this first example, 7/1 = 0. The source is located at a redshift of two, the 
loop at a redshift of one and the loop radius has been chosen to give an Einstein radius of Re/ R = 1 which we have 
shown to be a typical value. This choice corresponds to a loop radius of i? = 2.3 x 10^^ . The upper left window 
shows the mapping onto the source plane of a series of vertical lines spaced evenly at intervals of 0.02 R in the image 
plane. The upper right panel shows the critical curves, the images of infinite magnification, for this example. The 
lower left panel shows the unlcnscd images of a grid of source circles each defined by 15 points. Finally, the lower 
right panel shows the images observed for the loop lens along with the projection of the loop. (This field of sources 
is primarily illustrative; realistic galaxies are larger compared to realistic loops, and a subsequent set of figures will 
address this issue.) Images were found by the aforementioned method with a grid resolution of 200 x 200. Figure ^ 
shows the same two-loop now rotated through O2 = 0, but now we only show the lensed field and the critical curves. 
Figure ^ shows the same information as Fig. ^ but for a three-loop instead. The loop parameters in this case were 
chosen to match those used by Stebbins a = 0.5, sini/i — 0.5, 9i = 62 = 50 and -0 — 0. Finally, Fig. ^ shows the 
same three-loop, but now at a phase ^ — 120. 

In the next series of figures, we consider more realistic images that might be observed with a loop lens and a distant 
galaxy. Recall that we expect the typical loop radius to be TGiit/2TT. Taking F = 60 and locating the loop at a 
redshift of z = 1, one finds that the loop subtends an angle of 3.2 G^q arc sec, where G/^g = lO^G/^. We have already 
stated that G/ig 1 if strings are responsible for seeding structure formation, but we can more precisely constrain this 
value. The best data comes from the observation of the cosmic microwave background. If strings are responsible for 
the large scale temperature fiuctuations observed by COBE, then G/ig = 1.5 ± .5, as calculated by Bennett, Stebbins 
and Bouchet Below this range, strings cease to be interesting candidates for structure formation. Upper limits on 
the string tension have also been set by observing millisecond pulsars. The most recent data suggests that G/ig < 2 
We consider G/ig = 1.25 so that our loop diameter subtends an angle of 4 arc seconds. The angle subtended 
by the visible portion of a distant galaxy at redshift z '--^ 2, is roughly 1 arc second. In Figs. p[]lO|, we show several 
examples of the images formed with a circular source with radius 1 arc sec for each of our loop examples with R = 2 
acr sec. In each case the Einstein radius was fixed to be one corresponding to a source redshift of about two. Source 
locations were chosen in these examples primarily to show some of the more interesting features of the loop lenses. 



We will not explicitly consider loops that have a net velocity though our formalism applies directly to this case. We expect 
that the lensing by moving loops will look veiv similar to the lensing by stationary loops but that the loop will have a different 
effective shape due to the condition in eq. (p2[). 
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V. DISCUSSION 



The examples in the previous section may not represent an exhaustive sample of possible loop configurations, but 
they do show some of the generic features of loop lenses. For both the two-loop and the three-loop, we have considered 
configurations for which the projected loop lies near the Einstein radius and more compact configurations which reside 
well inside the Einstein radius. For the former case. Figs. |^, |^, ^ and ^, we notice that, like the perpendicular circular 
lens, we can have images which pass through the center of the loop and lie close to where the unlensed image would 
be. This is perhaps the distinguishing feature of loop lens images. In Figs. |^ and || we see examples of relatively 
undistorted images encircled by arcs and rings in a manner unlike what one would expect for a more homogeneous 
lens mass distribution. The more compact loops produce, unsurprisingly, less spectacular images. In Fig. ^ there are 
some examples with three images which may be distinguishable from ordinary lenses. Still, for all loops, one feature 
always differentiates them from ordinary lenses: namely that the lens itself is dark. 

Would the observation of lensed images lacking an observed lens in itself confirm the existence of cosmic string 
loops? Unfortunately, the answer is no; the existence of dark lenses in and of itself would not be conclusive since there 
are other possibilities. For example, it has been suggested that dark matter could form halos without a luminous 



component and these would result in dark lenses 1 18 19[| . Kandaswamy, Rees and Chitre |]19[ found that a single dark 
matter halo is not sufficiently dense to produce multiple imaging, but the alignment of two such halos could. Typical 
image separation was on the order of a few arc seconds, similar to the string loop case. However, these halos will have 
a distribution of matter consistent with coUisionless particles (i.e. an isothermal sphere), so their images should be 
distinguishable from string loops. A second possibility is that the net effect of many distant lenses which individually 
would not produce multiple images could add together to produce such an effect. This, however, has been shown to 
be statistically unlikely [gO[. Perhaps the simplest explanation for a dark lens is that it is not dark at all. A cluster 
of galaxies could produce two images which are separated on arc second scales, with a third located typically an arc 
minute away. The third image could easily go unnoticed creating the false impression of a dark lens pl| ]. But again, 
the images should be distinctly different than those of a string loop. The good news then is that string loop lenses 
are distinguishable from other dark lens possibilities, especially loops which are not compact with respect to their 
Einstein radius. Should a dark lens be confirmed, detailed observations of the image structure ought to determine 
the nature of the lens. 

Another possibility that we have not examined in detail in this paper but that is amenable to an identical analysis 
is the lensing due to wiggly long strings. The presence of wiggles on long strings means that the long string can 
probably be regarded as a sequence of small loops. In this case, we should observe a linear sequence of dark lenses 
with each lens having the characteristics of a string loop lens. Such a lens is unlikely to occur in the context of any 
other model. 



VI. CONCLUSION 



In this paper we have derived a method for calculating the deflection of light rays due to the gravitational field of 
an oscillating string loop. We have shown that this problem can be reduced to an effective static problem, greatly 
simplifying calculations. The formalism was then applied to the problem of gravitational lensing by cosmic strings, 
and using typical loop parameters, we have shown that a loop lens produces images on arc second scales, similar to 
galactic size objects. Specifically, we find that for Gfi ~ (1 — 2) 10~^ — values consistent with structure formation, 
microwave background anisotropics and millesecond pulsar timing limits — strings can produce images separated on 
arc second scales which would be observable by both ground based telescopes and the Hubble space telescope and 
would have features that are distinctly different from other dark lens candidates. This suggests that string loops can 
be definitively observed as gravitational lenses. Furthermore, the lensing due to long strings would appear like a linear 
sequence of lensings due to loops and these would be the unmistakable fingerprints of cosmic strings. 
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APPENDIX A: EFFECT OF A STRING LOOP ON THE ENERGY OF A PHOTON 



The temperature fluctuations induced by cosmic string loops has previously been investigated by Stebbins 
using the harmonic gauge as we have done here. Stebbins, however, used the real space retarded Greens function 
to find the change in 70, but, as we are solving the problem in Fourier space, it is useful to compare our results. 
Let us return then to eq. and replace ka with fc|| or equivalently kg. If we rewrite the product 1/0116*^11 as 
1/(1 - h)d/dt(eh'^^-f\\'>), we can again evaluate the fcii integral which now gives us 



da I dt I d^k, ^ '^^''^ 



dt 



(Al) 



The time integral can be evaluated by parts leaving 



da / d k 1 



-ik I f 



(A2) 



J t = to 



and we can again evaluate the integral as before. Finally, the temperature change caused by the loop is given by 



AT 



4G/Lt / da-f^-f" 



Fniy \ log(/2 ) 2F^lly fx • f± 



dt \ (I-/11) / (I-/11) (1-/||)2 fi 



(A3) 



The appearance of the logarithm in the first term in the above expression may be cause for concern as it appears to 
diverge for rays passing far from the loop. However, Stebbins |9| has proven, using conservation of energy that 



l^(i-/ii)j (1-/11) 



= 0, 



(A4) 



For rays far from the loop, i.e. xq ^ r, the logrithm can be expanded to leading order, log(/^) « log(a::Q) — 2xo • r/a;§, 
which falls like the inverse of the distance and does not diverge. 

For completeness, we sketch the proof of eq. (A4). Conservation of energy requires that d^T^^j = 0, so we can 
certainly write 



d-'x / (5(x • 7 ~ t)d''Tf,^{x, t) = 0. 



(A5) 



We are in essence performing the integral over the time slice as it appears in eq. ( |Al| ). Integrating the spatial 
components by parts, one can easily verify that eq. (Af;) is equivalent to 



d^'x I dtd{^-^~t)Y^T^,{^,t)^0, 



(A6) 



given that t = 7 • x is implied by the delta function. Using eq. (E2|) as a guide, we can write the stress energy as 



T, 



^ / daTf,^{a,t)S{-}<L- f{a,t)), 



(A7) 



and substituting into eq. (|A6|) gives us 



^l I da I d^'x I dt I dt' S{jc-j-t)Y^^[T^A^,t')S{^-f{a,t'))S{t-t')] = 0. 



(A8) 



We have added the extra delta function to enable us to see how to evaluate this integral. Let us replace it with a 
Fourier integral and evaluate the derivative. 



I da I d^x I dt I dt' I dfc 5(x-7-07%^(cr,t')'5(x-/((7,t'))^e*'''*'"*^ = 0. 



(A9) 



Now we can evaluate the integrals over t and x, leaving us with 
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da J dt' j dk 7-V(a,i')^e"=(''-^il(^'*'» = 0, 
which is equivalent to (rewrite the product of k and the exponential as a time derivative) 



fijdajdt'j dk -f" — 



d 



-5{ik{t' - f\\{a,t'))) = Q. 



This can be integrated by parts, finally giving us 



fi da 'f'^ 



T 



dt \l-f^^{a,t)J l-f\\ia,t) 



-0, 



(AlO) 



(All) 



(A12) 



where to is again the solution to /|| (a, to) ~ to- To recover eq. ( A4), one may contract the above with 7^ and recognize 
that 7^7"^P. = I'^l'^Ff,,- 

We should point out that both terms in eq. (A3) are necessary for calculating the microwave background anisotropics 
produced by string loops. In examples, however, we found that the contribution of the logarithmic term to the 
temperature fluctuation was only a few percent at most, and so we suspect that, in general, the logarithmic term that 
was overlooked in Ref . will not make a significant difference to the existing analyses . 
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FIG. 1. A schematic representation of a gravitational lensing system, 
is defined as the difference of unit vectors along SI and lO. 



The vectors rj, ^ need not be coplanar. The vector a 
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FIG. 2. The images resulting from a planar circular loop of radius one (solid circle) for a series of circular sources which 
spiral out from the origin. Hatching indicates images resulting from a particular source. The dashed line shows the Einstein 
radius selected for this example. 

FIG. 3. Lensing by a two-loop rotated through 9i = 0, 82 — 90 and with phase tp — 0. Each panel has dimensions 4ii x 4R 
where R is the size of the loop. The upper left panel shows the mapping of a grid of lines in the image plane onto the source 
plane. The upper right panel shows the critical curves in the image plane, i.e. the image points of infinite magnification. The 
lower left panel shows the unlensed images of a grid of circles each defined by 15 points. Finally, the lower right panel shows 
the lensed images along with the projection of the loop at time to- 

FIG. 4. The left panel shows the projected loop and images of a field of sources like those used in figure |^, but with a 
two-loop rotated through 9i — 62 — and phase ip = 0. The right panel shows the critical curve 



FIG. 5. Same as figure W, but with a three-loop rotated through Oi — 62 = 50 and phase tp = 



FIG. 6. Same as figure but with a three-loop rotated through 9i — 62 = 50 and phase tp = 120 

FIG. 7. The upper left panel shows the caustic of a two-loop rotated through 9i = 0, 62 = 90 and phase tp — 0. The other 
panels show the images as crossed points produced by a circular source whose unlensed image is given by the filled points. Also 
shown is the caustic. 



FIG. 8. Same as figure u\ but with two-loop lens rotated through 61 = 82 — and phase tp = 0. 



FIG. 9. Same as figure W but with a three-loop rotated through 9i — 62 = 50 and phase ip = 0. 



FIG. 10. Same as figure W but with a three-loop rotated through Oi — 62 = 50 and phase ?/) = 120. 
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